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Properly polarized and correctly timed femtosecond laser pulses have been demonstrated to create 
rotational states with a preferred sense of rotation. We show that due to conservation of angular 
momentum, collisional relaxation of these unidirectionally rotating molecules leads to the generation 
of macroscopic vortex flows in the gas. The emerging velocity field converges to a self-similar Taylor 
vortex, and repeated laser excitations cause a continuous stirring of the gas. 



PACS numbers: 33.80.-b, 37.10.Vz, 47.32.Ef 

Nowadays, femtosecond lasers are routinely used to 
manipulate and control molecular rotation and orien- 
tation in space [H (Sj. Particularly, diverse techniques 
have been developed to bring gas molecules to a rapidly 
spinning state with a preferred sense of rotation. These 
schemes include the optical centrifuge [3H5j and the 
molecular propeller [6l[7] methods, excitation by a chiral 
pulse train [8] and other techniques that are currently 
discussed [9j[T0]. In rarefied gas, coherent effects such as 
alignment revivals have been widely investigated. Most 
recently, the optical centrifuge method was applied to ex- 
cite molecules to rotational levels with angular momen- 
tum of hundreds of h in dense gas samples [5 , where col- 
lisions play an important role. In this paper we analyze 
the behavior of a dense gas of unidirectionally spinning 
molecules after many collisions have occurred. The sta- 
tistical mechanics and equilibration process of such a sys- 
tem are not trivial, as known for other ensembles in which 
angular momentum (AM) is a conserved quantity in ad- 
dition to energy [TTJ [12] . We show that due to the AM 
conservation, the collisions in such a gas transform the 
laser induced molecular rotation into macroscopic vortex 
flows. The lifetime of the generated vortex motion ex- 
ceeds the typical collision time by orders of magnitude, 
and the emerging velocity field eventually converges to a 
self-similar Taylor vortex [13 . The viscous decay of this 
whirl can be overcome by repeated laser excitations that 
produce a continuous stirring effect. 

Consider a gas sample excited by laser pulses causing 
the molecules to rotate, on average, unidirectionally. Our 
analysis of the subsequent gas dynamics was performed 
in two steps. First we used a Monte Carlo approach to 
simulate directly the kinetics of molecular collisions just 
after the excitation. Then, after the motion had taken 
a collective character, we used continuum hydrodynamic 
equations to investigate the gas dynamics and vortex flow 
evolution. 

We analyzed the initial stages of the transformation 
of the individual molecular rotation into collective gas 
motion with the help of the Direct Simulation Monte 
Carlo technique (DSMC [II]). This statistical method is 
a proven tool for computational molecular dynamics of 
transitional regimes [15 , although the standard DSMC 



does not always conserve AM [16l [17] . For this reason, 
we designed an advanced variant of the DSMC simulation 
in which angular momentum is strictly conserved. The 
simulation was run in axisymmetric geometry, and han- 
dled collisions using a simple 'hard dumbbell' model that 
treats each 'molecule' as a couple of rigidly connected 
hard spheres [18 . In recent experiments, laser excitation 
of unidirectional rotation induced an average axial AM 
in the range of (J^) ~ lOh - 400h [4, 7 . We first simu- 
lated a homogeneous sample of such unidirectionally ro- 
tating nitrogen molecules put at atmospheric conditions 
and confined to a smooth cylinder. In order to reduce 
the statistical fluctuations of the Monte Carlo simulation 
we averaged the results over multiple runs [19]. Figure 
[l] shows that within a few collisions (0.3 — 0.5 ns after 
excitation at our simulation conditions) the energy par- 
tition between the translational and rotational degrees 
of freedom approaches equilibrium. Despite the simplic- 
ity of our collision model, it provides the expected en- 
ergy partition ratio of 3:2 (translationirotation) typical 
for diatomic molecules. It takes about ten collision times 
for the molecules to completely lose memory of the pre- 
ferred direction of their rotation, supposedly even if the 
molecules remain in excited rotational states [5]. 

At longer times, the confined gas rotates rigidly (with 
tangential velocity linear with the radius) as seen in Fig. 
|2j This represents a steady state for the motion within 
a frictionless cylinder, eliminating shear stress. As ex- 
pected, the results show that the rotation speed is pro- 
portional to the total initially induced AM, predicting a 
steady macroscopic gas rotation frequency of ~ 10^ Hz 
for average initial unidirectional rotation of (J^) ~ lOOh 
at atmospheric conditions within a 10 /im diameter cap- 
illary. 

At the next stage we analyzed the formation of an un- 
confined vortex. For this we considered excitation by an 
axisymmetric laser beam with a Gaussian radial profile. 
Our study showed that with time, a circular flow develops 
around the beam's axis. We examined the spatial and 
temporal development of the circular flow structure by 
first running the DSMC simulation for the non-uniform 
initial rotational excitation. We sampled the results after 
thermalization, and used continuum gas-dynamic anal- 
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Figure 1. (Color online.) Mean translational energy, rota- 
tional energy and molecular axial rotational angular momen- 
tum evolution. Initially all the molecules were unidirect ion- 
ally rotating. The DSMC results show that the final transla- 
tionahrotational energy ratio reaches 3:2 (green:blue curves, 
right axis) within about 5 collision times. The average axial 
component of the molecular AM approaches zero at about 10 
collisions, showing that there is no longer a preferred rota- 
tion direction. Time axis is logarithmic in units of average 
collision time 10"^*^ s). 
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Figure 3. (Color online.) (a) Spatial dependence of the gas 
density and the tangential velocity of the unidirectional rota- 
tion at different times. The time (from left to right) is 8 to, 
15 to and 21 to after the simulation start (to is the typical vor- 
tex timescale, see Eq. ([3|). The emerging acoustic wave is 
clearly visible, and the lower density core can also be noticed. 
The color coding represents the tangential velocity value, also 
shown in (b) in the plane perpendicular to the laser beam at 
t — Sto. (c) Density radial profile at the same times. The 
spatial coordinates throughout are in units of wo, the radius 
of the beam's waist. 
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Figure 2. (Color online.) Tangential velocity profile of the 
gas, which exhibits rigid-body-like rotation. The data are 
averaged over multiple DSMC runs, covering a time range of 
t = 10~^ — 2 • 10~^ s. In this simulation, all the molecules 
initially rotate with AM of (Jz) ^ lOh. The bars indicate the 
statistical error of our procedure. The solid line is a linear fit 
with - 0.99 



ysis to unfold the subsequent gas behavior. The gas- 
dynamic equations we used relate to the mass, energy 
and momentum conservation in a compressible gas [20] . 
In order to better account for compressibility we added 
bulk viscosity [21 terms to the equations, and the ideal 
gas equation of state concluded the equation set. We 
used an axially symmetric finite-difference time-domain 
numerical scheme to solve the problem, with initial con- 
ditions based on the DSMC results. The full analysis of 
the numerical solution will be discussed elsewhere, here 
we report some of our observations. 

Since the laser injects energy to the molecules by ro- 
tating them, the gas at the beam's focus is heated follow- 
ing the laser excitation. The subsequent gas expansion 
creates a density crest, clearly visible on Fig. |3|i, that 
propagates outwards at about the speed of sound. 

The radial profile of the directed tangential velocity 
shows formation of a rigidly rotating core at the laser 
beam's focal spot. The profile exhibits a diffusion- like 
behavior in the radial direction. This phenomenon is at- 
tributed to the viscosity and is consistent with the notion 
of 'vorticity diffusion' [2T. After the emerging acoustic 
wave had left the core, naturally also the temperature 
spreads diffusively, and so do the other thermodynamic 
and flow variables (cf. the density plot in Fig. [s]^) . 

At that post-acoustic stage the radial profiles of the 
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Figure 4. (Color online.) Scaled tangential velocity ra- 
dial profile at different times. The tangential velocity Ue 



the radius r is scaled as 



is scaled as Uet^^^ ■ 7i ^ / 

rt~^'^^/[wot^^^^]. The plot shows a fast convergence to a self- 
similar Taylor vortex profile (shaded curve, red online) given 
by Eqs. The time t was offset by a vortex 'age' 

of 0.9 to to produce the best fit. 



velocity, density, pressure and temperature at different 
times seem to be self-similar. This can be illustrated 
by examination of the momentum conservation equation 
(tangential component): 

fdUe dUe ^ UrUe\ ^ ^ ( rj^\ 
\ ot or r J or \r or J 

where p, Uq^ Ur^l-i^ r and t are the density, tangential and 
radial velocity, viscosity, radius and time, respectively. 
The velocities (e.g., in Fig. ^p) are small compared to 
the typical velocity (Vb = ^/pwo ~ 100 ^/s, wq being the 
inducing laser beam's waist radius). We may thus ne- 
glect the second order velocity terms in Eq. ([T]). The 
relative density variations also shown in Fig. ^ are less 
than 10%, and by taking a nearly constant viscosity, Eq. 
([T]) reduces to a diffusion equation, with a 'diffusion co- 
efficient' proportional to the kinematic viscosity v = p 
[23] . Following that insight, we found self-similar fits to 
the radial profiles of the gas parameters (after the acous- 
tic wave has left the core). Specifically for the tangential 
velocity. Fig. [4] shows that the scaled calculated profiles 
converge with time to the Taylor vortex form [T3l [23] : 



Ue{r,t)=n 
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where the length scale wq is the radius of the waist of the 
laser beam, and the typical vortex lifetime is determined 
by to- The dimensionless parameter H depends on the 
average induced AM, and characterizes the number of 
turnovers the core makes before the vortex' viscous decay. 
The parameters to and H are defined as 



^0 = 



Wq 



n = 



STTnimol • ^ 



(3) 



where {Jz) and rrimoi are the average AM of the molecules 
excited by the laser and the mass of the individual 
molecules, respectively. 

It follows from Eq. ([2| that the rigid core of the vortex 
rotates at a frequency that drops with time as t~^, while 
the typical core size grows as t^/^. For a tightly focused 
laser beam (wq ~ 1 /im) that excites nitrogen molecules 
at atmospheric conditions to (Jz) ~ lOOfi, the vortex 
parameters are 1-L ~ 10~^, to ~ 10~^ s. In this case a 
substantial initial rotation frequency of ~ 10^ Hz may 
be achieved. 

Consider now the motion of gas parcels as the vortex 
develops and decays. We calculated the parcels' trajec- 
tories both numerically and analytically with the help of 
Eq. (§, mapping each point to its location after the 
vortex had died. The jolt imparted by the laser ex- 
citation results in a finite angular displacement of the 
gas parcel. By repeating the excitation again and again 
the angular displacement accumulates, rotating parcels 
around the beam's axis as a 'micro-torque-gun', produc- 
ing a motion similar to that of a ticking watch dial. At a 
pulse repetition rate of 1 kHz^ the time interval between 
pulses is long enough for diffusion to wipe the tempera- 
ture and density inhomogeneities and any radial displace- 
ment, even on millimeter-scale. Thus, in the micron-size 
waist of the laser beam only the azimuthal displacement 
grows from jolt to jolt. Figure [5] and the Supplementary 
Video [25] visualize the angular displacement after many 
iterations. This angular displacement pattern can be ob- 
served by following the motion of micron-sized particles 
suspended in the whirling gas. 

To conclude, we showed that laser pulses that induce 
unidirectional rotation of individual gas-molecules pro- 
duce a long-lasting macroscopic rotation of the gas. In a 
frictionless cylinder the gas reaches a steady rigid-body- 
like rotation. In free space a self-similar vortex forms 
and lasts significantly longer than both the pulse du- 
ration and the typical collision time. We also demon- 
strated that a repeated laser excitation can stir the gas 
in a pounding 'torque-gun' motion accumulating to large 
angular displacement. The 100 kHz rotational motion of 
the gas vortex may be detected with the help of the ro- 
tational Doppler effect [26H28] or by usage of structured 
light spectroscopy [29] . According to our estimations the 
predicted collective spiraling motion is also observable by 
monitoring the motion of particles suspended in the gas. 
Finally, we emphasize that the appearance of the gas vor- 
tex is a direct manifestation of angular momentum con- 
servation. In a sense, our laser-induced effect presents a 
table-top analogue to the formation of tropospheric cy- 
clones, in which global rotational structures appear due 
to the inverse cascade scale-up starting from the motion 
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Figure 5. Accumulated displacement deformation of a virtual 
rectangular mesh after 6,000 pulse sets, for 1-L ^ 10~^. The 
angular displacement is clearly visible around the center (cir- 
cle marker). The coordinates are in the plane perpendicular 
to the laser beam, in units of the laser beam's waist radius 
Wo. 



of turbulent eddies [30l |3T] . 
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